Stochastic Processes in Yellow and Red Pulsating Variables 
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ABSTRACT 

Random changes in pulsation period are well established in cool pulsating stars, in particular 
the red giant variables: Miras, semi-regulars of types A and B, and RV Tau variables. Such effects 
are also observed in a handful of Cepheids, the SX Phe variable XX Cyg, and, most recently, 
the red supergiant variable, BC Cyg, a type C semi-regular. The nature of such fluctuations 
is seemingly random over a few pulsation cycles of the stars, yet the regularity of the primary 
pulsation mechanism dominates over the long term. The degree of stochasticity is linked to the 
dimensions of the stars, the randomness parameter e appearing to correlate closely with mean 
stellar radius through the period P, with an average value of e/P = 0.0136±0.0005. The physical 
processes responsible for such fluctuations are uncertain, but presumably originate in temporal 
modifications of envelope convection in such stars. 

Subject headings: stars: oscillations — stars: variables: other — instabilities 



1. Introduction 

It is well known in variable star observations 
that it is almost impossible to predict the exact 
moments of light maximum in some late-type pul- 
sating variables, such as Miras and semi-regular 
variables. The cyclical light patterns displayed in 
such stars are reasonably well defined over long 
time intervals, but the regularity of their pulsation 
is typically marked by other effects best revealed 
through careful O-C analysis. A common effect 
is that of "random" fluctuations in pulsation pe- 
riod for a star from one cycle to anoth er. Eighty 
years ago lEddineton fe Plakidid (|l929h developed 



a novel method for establishing the importance of 
such random fluctuations for Mira variables, and 
the techni que has been revived frequently in re- 
cent years dPercy et alj|l997t IPercy fe Haldri 998: 



Percy fc Cohvaslll999l; IPercv et alJl2003L l2007h m 
order to establish its importance in other Mira 
variables as well as in other types of pulsating vari- 
ables. 

In essence the technique involves computing, 
without regard to sign, the average accumulated 
time delays (u(x)) between light maxima sep- 
arated by x cycles. If the time deviations in 
light maxima are dominated by random fluctua- 
tions in pulsation period, then as demonstrated 
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Fig. 1. — Schematic expectations for random pe- 
riod variations in a pulsating star with random- 
ness factor e = 2.5 and uncertainties of a — d .08 
in determining ti mes of light maximum, accor ding 
to predictions bv lEddington fe Plakidisl ( 1929t h 



by lEddington fc Plakidisl (|l929h the data for all 
available observed light maxima should display a 
trend described by: 

(u(x)) 2 = 2a 2 + xe 2 , 

where a represents the average uncertainty in days 
for established times of light maxima and e rep- 
resents the magnitude of any random fluctuations 
in period. Fig. Q] represents schematically such ex- 
pectations for a pulsating variable with a random- 
ness factor of e = 2.5 and uncertainties of a = 0.08 
day in established times of light maxima. 

2. Results for Different Pulsators 

Convincing evidence for random cycle-to-cycle 
variations in pulsation period for several differ- 
ent types of long-period variabl es has been gath- 
ered by Percy and c o llaborators ( Percy et al.lll997 ; 
Percy fc Halel Il998t iPercv fc Colivaj|l999l ) based 
upon the Eddington-Plakidis method. In contrast, 
mostly negative findings were obtained for shorter- 
period pulsators s uch as the (3 Cep star BW Vul 
(jPercv et al1l2003h and a sampl e of 6 Set and RR 
Lyr variables ( Percy et al.ll2007t ) , an exception be- 
ing the SX Phc variable XX Cyg, for which a weak 
signal of random fluctua tions in period was ob- 
tained (|Percv et al.ll2007l ). 

Similar searches have be en made using O-C 
data for Cepheid variables (jTurner fc Berdnikov 



200ll 12004 iBerdnikov et~aj]|2004 [20071. l2009alTbl: 
Abdel-Sabour Abdel-LatiJ 12004) . but generally 
with negative results in most cases (see results 
for j) Aql, S Cep, and T Ant in Fig. The 
long-period Cepheid SV Vul (Fig. O represents 
the case of a positive detection of random fluc- 
tuations, and also a case where one can observe 
individual light maxima. The situation is compli- 
cated for short-period Cepheids because separate 
light maxima are rarely observed and individual 
times of maximum often refer to data obtained 
over many adjacent cycles about the one cited. 
Since random fluctuations in period exist over 
several pulsation cycles, their effects on times of 
light maximum can easily be confused with other 
sources of scatter in the light curves. 
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Fig. 2. 
Cepheids. 



200 300 400 500 
x(cycles) 

Eddington-Plakidis tests for four 



Similar evidence for short time scale random 
fluctuat i ons in period can be seen in the results of 
Poleskil (|2008l ) for several LMC Cepheids. Offsets 
in the times of light maximum from cycle to cy- 
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cle are clearly present in Poleski's O-C diagrams, 
but the deviations are relatively small in compari- 
son with those observed over longer time intervals, 
where the effects of evolu tionary changes in m ean 
radius become dominant ( Turner et al. l2006ah . 

The M supergiant variable BC Cyg (M3 la) is 
an interesting example of random fluctuations in 
period for the SRC class of long-period variables, 
a small group of pulsating stars often overlooked 
by variable star observers. Times of light maxima 
for BC Cyg can generally be established to within 
a week or so (~ 0.01P), but the O-C variations 
(Fig. [3]) display scatter amounting to as much as 
a few hundred days. The best-fitting downwards 
sloped parabolic trend is otherwise consistent with 
the established period decrease in BC Cyg between 



1900 and 2000 ([Turner et aTll2006bh 
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O-C trends for the M supergiant BC 



A test of the O-C deviations for BC Cyg us- 
ing the Eddington-Plakidis technique is displayed 
in Fig. [H The deduced randomness parameter in 
this case of e = 12.88 is consistent with the rec- 
ognized dependence of the parameter on pulsation 
period P. In general, all existing results suggest 
th at the "e" parameter in creases with period P 



(jTurner fc Berdnikovll200lh . as shown in Fig. [5] A 
working relationship for the dependence is given 
by: 

e = -0.421(±0.146) + 0.015(±0.001)P . 
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Fig. 4. — An Eddington-Plakidis test for BC Cyg. 
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Random Period Fluctuations 
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Fig. 5. — The trend of increasing stochasticity e 
with pulsation period. 

3. The e/P Parameter for Stochastic 
Changes 



It was noted by lEddington fc Plakidisl (|l929l ) 
that the observed random fluctuations in period 
for o Cet (Mira) and \ Cyg amounted to 1.35- 
1.39% of the pulsation period, a parameter that in 
turn depends directly on stellar radius. In other 
words, a better parameter for describing stochas- 
tic processes in pulsating stars should be the ra- 
tio e/P, which must be independent of radius if 
the variables obey a period-radius relation. Re- 
sults to date for all Eddington-Plakidis analyses 
of pulsating stars (shown in Fig. [6]) confirm that 
assumption. The parameter e/P is indeed rela- 



3 



tively independent of pulsation period, i.e., inde- 
pendent of stellar radius, and has a mean value of 
0.0136 ± 0.0005 (±0.0 69 s.d .), identical to what 
Eddington fc Plakidisl (| 19291) concluded 80 years 
ago. 
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Fig. 6. — The ratio e/P appears to be independent 
of pulsation period. 



4. Conclusions 

In stars such as SV Vul (Fig. 0) and BC Cyg 
(Fig. [4]) the computed values of (u(x)) 2 initially 
increase directly in proportion to increasing differ- 
ences in cycle count x, cis p redicted for a stochastic 
process ( Eddington fc Plakidisl 19291 ). At larger 
cycle differences, however, the trend is reversed as 
the dominant pulsation in such stars reimposes its 
regularity in the observed times of light maximum. 
The observed values of (u(x)) therefore become 
very small for large cycle differences. Similar char- 
acteristics are observed in m any different pulsating 
stars, and was also noted bv lEddington fc Plakidisl 
(<1929l) . Stochastic fluctuations in period appear to 
be a common feature of nearly all pulsating stars, 
but they are invariably dominated by the regular 
pulsation in such stars. The physical processes re- 
sponsible for such fluctuations are uncertain, but 
presumably they originate in temporal modifica- 
tions of envelope convection in such stars. 
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